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Анотація. У статті розглядається задача оптимізації топології виробів у аддитивному виробництві 
завдяки оптимальному розміщенню кругових отворів. Задачу зведено до пакування кількох кругів змінних 
радіусів, заданих у | межах, встановлених стандартами  3Ю)-друку. Запропоновано  двокритеріальне 
формулювання, в якому враховуються коефіцієнт пакування та максимальне механічне напруження виробів. 
Для пошуку компромісного розв'язку задачі застосовано метод головного критерію. Розроблено новий підхід, 
який грунтується на модифікованому методі Аполонієвих пакувань кругів та нелінійній оптимізації. Наведено 
числові приклади та графічну ілюстрацію результатів. 


Ключові слова: аддитивне виробництво, оптимізація топології, двокритеріальна задача, метод головного 
критерію, пакування кругів, нелінійна оптимізація. 


Вступ оптимальної геометричної форми та 

Однією з найважливіших тем топології розробленого виробу, що 
сучасного машинобудування є зниження забезпечує мінімальну вагу за заданої 
ваги | за | збереження специфічних міцності. З цією метою для пошуку 
характеристик конструкцій, що найліпших параметрів конструкції, що 
проєктуються. Такі задачі з суперечливими задовольняють технологічні Й механічні 
критеріями пов'язані з прийняттям рішень обмеження, | використовують методи 
у сучасному виробничому процесі. оптимізації топології (1). Під час 
Ключовою метою задач | є пошук оптимізації топології конструкції рівень 
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механічних напружень у певній частині 
конструкції може бути використаний як 
індикатор неефективного використання 
матеріалу. У ідеалі рівень напруження має 
бути | рівномірним | і | близьким до 
граничного. 

Застосування методів оптимізації 
топології у машинобудуванні - це нова 
розробка | в | процедурі  проєктування. 
Найбільший імпульс у своєму розвитку ці 
методи отримали, коли стало можливим 
використання аддитивних технологій у 
процесі виробництва |3| замість класичних 
методів. Аддитивні технології дозволяють 
розширити асортимент конструкцій. 


У останні два десятиліття 
оптимізація топології стала активним 
полем для наукових досліджень. Це 


призвело до використання 
багатодисциплінарного підходу у пошуках 
вирішення проблем, в якому одночасно 
використовують методи механіки твердого 
тіла, термодинаміки та біології |4|. 

У роботі (|2| розглянуто 
методи оптимізації топології. 

На сьогодні можна виділити такі 
основні методи оптимізації топології: 
ЗІМР (метод твердого  ізотропного 
матеріалу З пеналізацією), Е5О 
(еволюційна структурна оптимізація), 
І еуе!-5еї (метод встановлення рівня) та їх 
різні комбінації. Одним з найбільш 
ефективних застосувань цих методів є 
оптимізація топології суцільних структур, 
тобто пошук найліпшого місця 
розташування та геометрії для отворів 
(порожнин) у зоні моделювання. 

У цій статті г ми | пропонуємо 
двокритеріальне формулювання для 
оптимізації топології на основі кругових 
отворів, яке враховує як коефіцієнт 
пакування, так і максимальне механічне 
напруження виробів. 

Задачі пакування кругів (СРР) мають 
широкий спектр | промислових та 
інженерних застосувань, | наприклад, 
проєктування розміщення об'єктів, задачі 
складування | в скляній, металевій, 
паперовій, текстильній та деревообробній 
промисловості тощо. 

Як правило, задачі | практичної 
оптимізації мають кілька суперечливих 


сучасні 
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критеріїв. | Такі | задачі | називаються 
багатокритеріальними. Методи оптимізації 
для задач багатокритеріальної оптимізації 
та її застосування розглядають у |, 6). 

У цій в роботі ми | пропонуємо 
двокритеріальну формулу для оптимізації 
топології на основі кругових отворів, яка 
враховує як коефіцієнт пакування, так 1 
максимальне механічне напруження на 
вироби. Запропоновано новий підхід до 
пакування кругів на основі 
модифікованоого скінченого Аполонієвого 
пакування кругів (АПК) та оптимізаційних 
гомотетичних перетворень кругів. Обидва 
прийоми дозволяють вигідно розмістити 


круги в багатокутних частинах, 
мінімізуючи матеріальні витрати. Для 
знаходження | компромісного | рішення 


застосовують метод головного критерію 
А 

Відповідно до АПК (див., наприклад, 
І31), починаючи з трьох взаємно дотичних 
кругів, наступні круги додають так, щоб 
вони були дотичними до трьох кругів, які 
утворюють чотири взаємно дотичних 
круги (або дотичні до межі контейнера). 
Ми  модифікуємо | АКІЇ, враховуючи 
верхню межу радіусів кругів та мінімально 
допустиму відстань між кругами. Деякі 
круги можна переміщувати всередині 
контейнера й розміщувати в місцях, які 
дають змогу збільшити радіуси наступних 
кругів. 

Методи аналізу напруженого стану в 
частинах З різною геометрією 
обговорюються в |9, 10). 

З кожним роком кількість робіт, що 
стосуються СРР, дедалі зростає. У |(11| 
подано огляд вибраних робіт, присвячених 
методам пакування та застосувань СРР. 
Зазначимо деякі останні статті (12-14. 
Потужним інструментом вирішення СРР з 


нерівними | кругами || є врахування 
додаткових змінних метричних 
характеристик кругів та/або контейнерів 
(див., наприклад,  |13-16)). | Найкращі 
результати | для | набору | еталонних 


прикладів СРР постійно оновлюються на 
відомому вебсайті (17). 

У роботі (10) запропоновано модель 
пакування кругів та швидкий евристичний 
алгоритм для оптимізації геометрії деталей 
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з урахуванням обмежень РМІ,5 (пряме 
лазерне спікання металів). Особливістю 
моделі є те, що радіуси кругів визначено 
заздалегідь. 

У П || представлено модель та підхід 
до чисельного рішення щодо пакування 
яйцеподібних об'єктів довільного розміру 


та орієнтації в оптимізовані опуклі 
контейнери. Площу контейнера 
мінімізують. Стратегія рішення 
грунтується | на аналізі | множників 


Лагранжа та нелінійній оптимізації. У 
рамках універсальної моделі круги та 
еліпси вважають окремими випадками 
яйцеклітини. 

Метод головного критерію |7| часто 
використовують для вирішення задач 
багатокритеріальної оптимізації. У цьому 
методі оптимізацію проводять лише за 
одним критерієм, водночас на решту 
критеріїв накладають обмеження. У |З| 
досліджують оптимальність рішень 
Парето, отриманих методом головного 
критерію. 

Двокритеріальну функцію для 
пакування об'єктів у модулі космічного 
корабля запропоновано в |18|. Розміри 
модуля зведено до мінімуму, а також 
враховують інші критерії, зокрема бажане 
сусідство між предметами та витрати на 
пакування. 

У роботі 
генетичний алгоритм 
двокритеріальної топології: мінімізують 
масу та  максимізують ефективність 
жорсткості на вигин та кручення. 


19) використовують 


оптимізації 


Постановка задачі 
Нехай | Р, є  багатокутником з 
Уді З (Хді» ді, Є Їо 2 Ц,2..55ої. 
Визначимо незв'язну багатокутну область 
Р-ЦОВ, 
ІєЇ,, 
РІРоФІєЇєІТ, «-ЩОоотьї, 


вершинами 


Іє І 


РІС Ро, ря 


де 

Руз( (5 у) є В" Фу у) 20, 

т 1,2, МІ) 
є опуклими багатокутниками з вершинами 
ур З (р Ук), 1ЄД з (1,2... 
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Фі (», у) - алі Я ФС -0 , 


нормальними рівняннями ребер Р, ІєЇ,. 


РЕЯ за з 


Нехай також визначено набір кругів 
СС) є (у є ВУ (хх 


НузУр єп) 


із змінними радіусами г 


з Та вектори 


трансляції уд З Ср» г)» 


чеїуніаа Мч оковр'єр, г, 


г" є нижнім та верхнім дозволеними 


ач 


значеннями ДЛЯ г рочз 


б максимально 


можливий радіус круга, в який можна 
вписати Р, ІєЇ,. 

Умови пакування кругів С, (у, Г,), 
де Їх, в областьР визначають наступним 
чином. 

Розміщення С, (удг,) У Р 


сиру є Рв іні Є, (у, г)Ї Р'-2, ЦП) 


4 Їх, 
де Р'- У" ЛіпР. 

Мінімально допустима відстань між 
Сара? та Со (ру з 


аізС(уг Є (у от, ))2рР, (дьє)є9,, (2) 
ге зро 


де 
фі5С, (уд ту С, (Уго) З 


р(а, Б)? 


ш тіп 
аєСд (ду), БєС, (уозу? 


р(а,Б) - це евклідова відстань між точками 
а,реє уз Й 
9; «Газ: СУ пос В, 
СуУзтєЄб КР) 


Задачу двокритеріальної оптимізації 
топології кругів в багатокутну область 
формулюють у такий спосіб. 

Розмістити круги С, (У,»Г,)» 4ЄЇу, 


у області Р, забезпечуючи розміщення 2, 
кругів у Ріх  ФЕба 
а м епеМ та обмеження на відстані 


р 
(2), з | максимальним | коефіцієнтом 
пакування та мінімальним максимальним 


механічним напруженням. 


багатокутник 
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Умови пакування (1), (2) аналітично 
описані за допомогою методу Ф-функції 
120), що дає змогу подати математичну 
модель як таку двокритеріальну задачу: 


тах к(о), -с(0)) 


3) 


за умови о є //, 


де зе(у,г) є 


ЗМІННИХ; у - (у, У2аазУ,) Є ВектОром ЗМІННИХ 


вектором 
параметрів розміщення; г -(п, аг) з 
вектор змінних радіусів кругів; функція 


к(ф)-л еа ї 


4єї, 


(4 


- загальна сума площ кругів; с(с) Є 
неявною функцією, визначає 
максимальне механічне напруження 
залежно від вектора у координат центрів 


кругів 1 вектора г радіусів кругів; 


яка 


У зіоє во": фик("рУєтрід-р20, 


(4,2)є 0), І є І, Фо (уд) 20, 4єЇ,» 


й ен 20,4єЇ,, з ки 20,4є1,) 
є областю допустимих рішень; 
2 
4в (УрУроТрїзд З ро хо 
2 2 ' 
«орто ЗК ер) 

є | псевдонормалізованою 
кругів Са та С, , (42)є 0); 


Фі (уд) ві, 


Ф-функцією 


піп 


Ф (у,,г)-тахі 
47 474 тіа,2,..М 


їєї, 
є Ф-функцією круга С,, 4є1Ї,, і множини 
Р зВ'МіпіР. 
Нерівність Фих (оз УдзТаь То 20 
забезпечує умову (2) для кругів С, та С,, 
(48)є9,, 1 нерівність  Ф,(у,,г,)20 - 
умову (1) для розміщення круга С, у 


багатокутнику Р. 


Функції к(с) та - (05) Є 
конкурентними. Отже, наша мета - знайти 
компромісне рішення  двокритеріальної 


задачі (3) - (5). 
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Підхід до вирішення задачі 

Для вирішення задачі (3) - (5) ми 
використовуємо стратегію мультистарту 
121). 

Стратегія вирішення. 

За заданого значення максимального 
механічного напруження ми вирішуємо 
задачу за методом головного критерію |71, 
в якому задачу (3) - (5) зводять до 
однокритеріальної 

пазах к(с5) (6) 
за умови оєУ/", 
У" «(ее |-с(о)чнє20). 

Основна мета - максимізувати суму 
квадратів радіусів кругів, які потрібно 
розмістити. | Неперервне відстеження 
механічних | напружень є складною 
задачею. Важливе значення має граничне 
значення є механічних напружень. 

Тимчасово послабимо нерівність 

-с(бф )чє20 (7) 
5) 

Тоді можна апостеріорно оцінити 
максимальне механічне напруження для 
отриманої топології області 


ву ман | Є, (9') та перевірити умову (7) у 


точці локального максимуму б). 

Якщо умова (7) є справедливою, 
тобто -с(ф)-є20, тоді існує допустима 
точка задачі (6) - (7). У противному разі 
обчислимо іншу точку локального 
максимуму с і повторюємо процедуру, 
доки не буде виконано (7). 

Задачу (6) - (7) зведемо до / підзадач, 
окремо для кожного багатокутника Р., 


І є Т, - Наша мультистартова стратегія має 
такі основні етапи. 

Етап 1. Використовуючи алгоритм 1 
на основі пакування рівних кругів (141, 
визначають нижню межу 2, кількості 2, 
кругів, | які | можна розмістити | у 
відповідному багатокутнику Р, ІєїЇ, З 
максимальним Будують 
точку и. 


коефіцієнтом. 
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Етап 2. За алгоритмом 2 на основі 
АПК обчислюють допустимі точки для 
задачі (6), починаючи з точки и, (8). 
Будують точку му) . 

Етап 3. Використовуючи ІРОРТ |22,|, 


обчислюють точки локальних максимумів 
задачі (6) - (7), починаючи з точки 


ф'з(и Мн М) 


Етап 4. Вибирають найкращу точку 
локального максимуму, що задовольняє 
нерівність (7), знайдену на етапі 3. 

Детально розглянемо нашу стратегію 
рішення. 


Алгоритм 1. Пошук нижньої межі 
кількості кругів 

Результати пакування 
залежать від кількості 
початкових точок. 


сильно 
кругів та 
Ми використовуємо 
ідею АПК. За даними АПК, круги 
розміщують таким чином, щоб вони 
торкалися трьох інших кругів |8|. Для 
роботи | АПК будують початкову 
конфігурацію дотичних кругів. Для цього 
реалізують пакування рівних кругів з 
радіусом г" (якщо це є можливим). Задачу 
класифікують як ПРР (проблема пакування 
ідентичних кругів) 0 |23|. | Алгоритм 
грунтується на моделях, описаних у | 14). 
Ми розглядаємо пакування кругів у 
багатокутнику Р, Іє І, . Нехай розміщено 


пре бр і; кругів в багатокутники Р), 
е2ї-Ї и 


дапі і. За умовою, г, обмежене зверху 


Тепер розмістимо | С 


г", що може бути загалом меншим за 


максимально можливий радіуст,, 1єЇ,, 


круга С, (Уд»ті ), вписаного в РІ. За таких 


умов центру, круга С,(у,,г') може 


вільно рухатися в багатокутній множині 
Р,с Р, (рис. 1), задає 


координати центрів С (у, ) У Рі. 


яка можливі 
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Рис. 1. Розміщення першого круга в Ві 


Обирають випадкове положення 
центра круга, а потім виконують 
покрокову процедуру послідовного 


додавання наступних кругів, відповідно до 
алгоритму, описаного в | 14, 24). Якщо круг 
з радіусом г" не можна розмістити в Р), 
переходять до етапу 2 (АПК), починаючи з 
круга з максимально можливим радіусом 

ж . 
п, ІєЇ,. Якщо для круга є кілька 


можливих місць, вибирають місце, для 

якого круг є дотичним до максимальної 

кількості ребер багатокутника. 
Розв'язують задачі 


тах гук» ке 0,12, 
(8) 
шо а у 


(ю 
де и! 


(Уві Є с В), якщо к-0 та 
й» - (Ха» Уазз ню Уачію Гачк є У? о во 
для пд КЛ 

М а шо ЄВРО пу у б Уд) 20, 

т -1,2,...М| 

реїоз 


9 


(хи | Яні - ро - зані - 
(ак ур) 2-0), 
ііз12,.К іс, ЕєіЗ,..К), 
(хижі "Хачнк у нія 7 Удчк )? - 

2 
"Ся пи за) 20, 


іа12 ак, ВЕЇ2 КІ 
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Обчислюють локальний максимум 
(9) - (9). 


продовжують розв'язання задачі (3) - (9) 
для наступного значення К, враховуючи 


. ж 
задачі Якщо гир - г, то 


ж 


сою ж ж ж 
и з (ха Ура жрносі» Урі» Хоніо Урі 0)» де 


ХамоУдчк  Обирають випадковим чином 
ж - ж й У М 
(ці є МО). Якщогук «го, ітераційний 
процес зупиняють. 
Точка 
Уге(у,п)е (ха, Удь оз Х дної» У дикі» 
Хак УаньГ,..Г , Гаанкі 


Ка 
1-1 
- це результат роботи алгоритму. Тоді 
каб, круг Січ Торкається принаймні 
1 


трьох кругів (або межі Р,) 1 він є першим 
колом, відповідно до АПК (рис. 2). 


Рис. 2. Побудова конфігурації, відповідно до АПК 


Алгоритм 2. Побудова допустимих 
початкових точок 


Ж 
Якщо ги 


максимумом задачі (6) - (7), то круги 


є строгим локальним 


Са» Сабо Є жорстко закріпленими, 


утворюючи «заклинене» пакування (25. 
Відповідно до АПК, розміщують наступні 


круги ка Сіно" Санні» починаючи 
з радіуса Го нрчі ДО каз» 7 КІЛЬКІСТЬ 


додаткових кругів після роботи АПК 
(рис. 3). Таким чином, отримують точку 


ми я (Ха » Уд оо» Ханмі-1, УдМ-1» 
з ж ж ж 
паці скін НО 
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Загальна кількість 


кругів, 
- ж 
озміщених УР, є Жак ч-к,, Значення 
УГГІ І 1 А5Р, 


радіусів належать сегменту |" ,г |. 


Локальна оптимізація 
Після побудови початкового 
розміщення для кожного багатокутника 


Р, ІєЇ,, загальна кількість кругів - 


ве а, М , 


Переходимо до задачі (6) - (7). 
Обчислюємо точку локального максимуму, 
а потім будуємо початкову точку 


о зи мом, )- Для розв'язання задачі 
нелінійного програмування  використо- 


вують розв'язувач ІРОРТ |23| разом зі 
стратегією декомпозиції (13). 


ЧК Ндср 


Рис. 3. Розміщення круга, відповідно до АПК 


Обчислюють декілька точок 
локального максимуму та 
впорядковуються у порядку зменшення 
к(о) (див. задачі (6) - (7)). Якщо немає 
точки, що задовольняє граничне значення 
с, (див. (7)), обчислюють нові точки 
локального з максимуму або обирають 
компромісне значення с, - Дс . Розрахунок 
значень механічних напружень у цій 


роботі не розглядається. Ми посилаємося 
на (10). 


Як показують обчислення, точки 
локальних | максимумів близькі | або 
збігаються З вихідними точками, 


побудованими відповідно до АПК. 
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Числові результати 

Ми розглядаємо контрольний 
приклад, представлений у (10). Досліджено 
залежність кількості розміщених кругів та 
щільності пакування від мінімально 
допустимої відстані між кругами. 

Приклад І.Р, задано 11 вершинами 


Уді 7 (0,0), У» 7 (0,9), Убз - (18,9), 
Ус 2 (0,28), Удз - (0,33), Удз - (0,40), 
Удв - (0,40), Удт - (60,40), Удв 7 (76,34), 
о З (98,11), му є (100,6), 


дл 71000). РА иа де Їь он), 


вершинами у|; 7 (22,31), 
УІ2 - (35,21), УМ3 г (6,8); 
Узі 7 (40, 25), уз» - (54,9), ущз - (12,7); 


РьфебЇзі- Є 


Уз1- (44,30), У32 - (68,12), У33 - (59,8), 
Узд - (42,28); уді 7 (42,36), Удз з (63,36), 
уд з (89,30), || узд є (73,13); | уві « (74,37), 


Уго - (95,37), узз - (91,31), узд - (74,35). 
Мінімальна дозволена відстань між 
кругами становить ро0; г «0.5; г' «5; 


кругів, 
п.91 


Загальна кількість 
розміщених уР Є 
(М 215,0 є 23,13 214, е31,3., 28). 
Значення критерію у вихідній точці (АПК) 
к(об) -361.3287, точці 


локального максимуму -- к -362.0946. 
Ілюстрацію розміщених кругів наведено на 
рис. 4. Цей приклад має значення для 
максимізації коефіцієнта пакування і немає 
сенсу через нескінчену механічну напругу 
зараб. 


становить а в 


Приклад 2. Див. приклад І. р- 0.5. 
Значення критерію у вихідній точці - це 
к(од) -317.9666, а у точці локального 
максимуму - к -319.3224. 


показано нарис 5: п-обб 
ме, М з 19, 24 210, 2 221, 2 25). 


Ілюстрацію 


Рис. 4. Розміщення кругів, 
відповідно до прикладу І 
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Рис. 5. Розміщення кругів, 

відповідно до прикладу 2 

Приклад 3. Див. приклад І. р- 0.75. 

к(од)-301.0078, к -301.7119. Ілюстрацію 

показано на рис. б: п-50 (М -8, 2» -12, 
23 28, М 17, Аз 5). 

Приклад 4. Див. приклад І. 


раї. 

к(оо) -288.9095, к -289.7335. Ілюстрацію 
показано на рис. 7:п-47 (М -8, 2» - 12, 
ха 28, Хд н14, Зо 25). 

Для всіх прикладів час виконання 
становить 30 секунд. 

Як показують результати обчислень, 
зі збільшенням мінімально допустимої 
відстані між колами внесок великих кіл у 
загальну площу кола збільшується, 
оскільки він пропорційний  квадратам 
радіусів кола. Очевидно, це призводить до 
зменшення кількості кіл з малими 
радіусами та збільшення радіусів інших 
кіл. 


У зв'язку з цим, вибір р має значення. 
Вибирається | компромісне значення: 
занадто низьке значення може призвести 
до низької якості друку через технічні 
особливості процесу, а надто велике 
спричиняє додаткову витрату матеріалу. 
Крім того, другорядна метас(а) на нього 


може впливати зміна значень р, г таг"'. 


Рис. 6. Розміщення кругів, 
відповідно до прикладу 3 


Рис. 7. Розміщення кругів, 
відповідно до прикладу 4 
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Висновки 

Оптимізація топології Є ключовим 
моментом в аддитивному виробництві. 
Моделі макетів із круговими отворами 
допомагають скоротити матеріальні 
витрати, відповідаючи вимогам міцності. 


У цій роботі досліджено 
двокритеріальну | модель оптимізації. 
Запропоновано | швидкий ефективний 
алгоритм пакування, який грунтується на 
Аполонієвому пакуванні кругів, 
нелінійному програмуванні та методі 


головного критерію. Підхід дає змогу 
оцінити кількість необхідних кругових 
перфорацій та знайти розв'язок задачі, 
максимізуючи коефіцієнт пакування та 
дотримуючись граничного | значення 
механічних напружень. 

Подальші дослідження можуть бути 
спрямовані на одночасний аналіз площі, 
охопленої | кругами та механічних 
напружень під час процесу оптимізації з 
урахуванням умов балансування 


21,26, 271. 
Робота | частково | підтримується 
«Програмою державних пріоритетних 


наукових досліджень та технологічних 
(експериментальних) розробок Відділення 
фізико-технічних проблем енергетики 
Національної академії наук України» та 
Національним науково-дослідним фондом 


України (02.2020/167). 
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